Topological Quenching of Spin Tunneling in Mni2-acetate 

Molecules 



Chang- Soo Park* and Anupam Garg 
Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208 

(February 1, 2008) 



Abstract 

We investigate the topological quenching of spin tunneling in Mni2-acetate 
molecules with an applied magnetic field along the hard axis. The model 
Hamiltonian describing this system has a fourth-order term due to the tetrago- 
nal anisotroy. We treat this model using the discrete phase integral formalism 
in which the Schrodinger equation corresponding to the Hamiltonian becomes 
a nine-term recursion relation. We solve this recursion relation and find that 
the tunnel splitting is quenched at certain values of the applied field. We also 
present a qualitative treatment based on the instanton approach. 
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I. INTRODUCTION 

Quantum tunneling of a spin or spin-like degree of freedom has been discussed for over a 
decade now |]l|, but unambiguous evidence for its existence has only come recently [Q] from 
studies on the magnetic molecule [Fe802(OH)i2(tacn)6]^"*' (Fcg for short). This molecule has 
a total spin 5 = 10, is biaxially symmetric, and can be modeled by the spin Hamiltonian 

n = hSl + k2Sl - g^iB^ ■ H, (1.1) 

with ki > k2 > 0. In the first approximation, it has two degenerate ground states, approxi- 
mately given by Sz = ±10, which are separated by an energy barrier in the xy-plane. The 
question of interest is to understand how these states are admixed by quantum tunneling. 

Direct numerical diagonalization of Eq. (|1.1|) using the experimentally determined values 
ki ~ 0.33 K, k2 ~ 0.22 K H-g reveals that the tunnel splitting A is ~ 10"^ K, which 
is too small to be observed directly. Wernsdorfer and Sessoli overcome this difficulty 
by applying a small amplitude ac magnetic field along the z direction, which causes the 
Sz = ±10 levels to cross one another. Transitions between these levels are now possible 
via the Landau-Zener-Stiickelberg (LZS) process 0, and the underlying tunneling matrix 
element A can be deduced from a measurement of the incoherent LZS relaxation rate for 
the total magnetization. The key experimental fact that supports this interpretation of 
the relaxation (which could after all be due to a classical activation process a priori), is a 
systematic and remarkable oscillatory dependence of the inferred splitting A on the strength 
of the magnetic field H when this field is applied along the hard direction x. This phenomon 
was predicted some time ago [0 based on an instanton approach. Briefly, when H || x, there 
are two symmetry related instanton paths that wind around x in opposite directions, and 
together form a closed loop on the (complexified) unit sphere. The actions for the instantons 
are complex, and differ by a real valued Berry phase given by S times the area of the loop, 
giving rise to interference ||^. Now, however, the Berry phase is not fixed at 271 S, but may 
be continuously varied by varying H. As a function of H, therefore, the tunnel splitting 
oscillates, and is completely quenched at values of H where the Berry phase is an odd integer 
times vr. 

In this paper, we wish to study topological quenching of tunneling in a second magnetic 
molecule, Mni2-acetate (or Mni2-ac for short), which has also been the subject of several 
experimental studies ||9|-[r3|. The reason for our interest in Mni2-ac is that, in contrast to 



Fcg, it has tetragonal symmetry. The spin Hamiltonian of Mni2-ac can be written as 

n = -ASl-BSt + C{Sl + St), (1.2) 

where A 3> -B ^ C > 0. The easy axis z now has four fold symmetry, the hard axes are ±x 
and ±y, and the medium axes are the lines y = ±x in the xy-plane. Here, the symmetry of 
a pair of instanton paths is preserved when a magnetic field is applied along one of the four 
hard axes. Thus, the quenched spin tunneling phenomenon is also anticipated in Mni2-ac. 
Whether it can be observed or not depends on how strong the environmental decoherence 
is, and is a question that we shall not investigate here. 

To investigate the topological quenching of the spin tunneling in Mni2-ac molecules we 
use the discrete phase integral (DPI) (or, discrete WKB) method |I^JT5[] . The DPI method 



has been applied to spin tunneling problems in a recent series of works by one of us (AG) 
1^19|. This method is semiclassical just as the instanton approach is, but it is easier to 



use for the study of the splittings of higher pairs of levels. It is particularly well suited to 
study tunneling when H is not along the hard axis. As in Feg, we anticipate that A will be 



quenched at a number of such field values [T^JS^ , which correspond to non-trivial diabolical 



points |]21| , P2| in the magnetic field space. In this paper, however, we shall not consider such 
general orientations of the field as the calculations for H || x are already quite complex. 

In the following section we present the DPI formalism for the present model. Unlike 
the previously studied model for Feg, the Schrodinger equation corresponding to the spin 
Hamiltonian for Mni2-ac becomes a nine-term recursion relation because of the fourth-order 
term. We give a systematic analysis for this recursion relation. We then calculate the tunnel 
splittings as a function of the applied field for the first few energy levels. The results will 
be compared with those obtained by numerical diagonalization of the spin Hamiltoian. In 
Sec. Ill we give a qualitative discussion of spin tunneling based on the instanton approach. 
This gives a good physical picture of the quenched spin tunneling and explains some of the 
interesting features found in the DPI results. A summary of the results in Sec. IV concludes 
the paper. 

II. DPI CALCULATION OF TUNNEL SPLITTINGS 



We consider the spin Hamiltonian in Eq. (|1.2| ) with magnetic field applied along the x 
axis. For convenience we divide the Hamiltonian by A to work with dimensionless quantities. 
With this choice we can write 

n = -si - x^st + USX + St] - sh,s,, (2.1) 

where Ai = B/A, A2 = C/A, K = H^/SHc {He = A/q^ib), and 

S = S + ]^. (2.2) 



Here, /i^ is the Bohr magneton, g = 2, and S is the spin. Following Ref. [T^, A/ks = 0.556K, 
B/kB = 1.1 X lO-^K, C/kB = 3 X IQ-^K, so that Ai = 1.98 x 10-^ A2 = 5.4 x 10-^ and 
He = 0.414T. Let |n) = \6, 0) be the spin coherent state with maximal spin projection along 
the direction fi, with spherical coordinates 6 and 0. We introduce the classical energy 

ncie,(i)) = {h\n\h) 

= -S"^ cos^ e-XiS^ cos^ e + 2X2S^ sin^ 6 cos 40 - SSK sin 6 cos 0. (2.3) 

When hx = 0, Tic has minima at ^ = 0, ^ = vr. As h^ is increased, these minima move 
toward 6 = 7r/2, = 0, lying in the a:2;-plane. At a certain critical field, hxco, these minima 
will merge with each other, giving rise to a double zero of d7ic{0, = 0)/d6 aX 6 = 7t/2. By 
using this condition, we can show that 

2 S* 
/ixco = ^(l + 4A252). (2.4) 



With the experimental numbers given above, h^co = 1.946 [|23 |. 



A. Recursion Relation 

The DPI formalism can be started with Schrodinger equation in the Sz representa- 
tion. Introducing TClip) = E\ip), Sz\m) = m\m), {m\ip) = Cm, and {m\H\m') = tm,m', 
the Schrodinger equation corresponding to the Hamiltonian (2.1) can be expressed as 

m+4 
n=m— 4 

This is a nine-term recursion relation with diagonal terms tm.m from S"^ and S^, and off- 
diagonal terms tm,m±i, tm,m±,4 which are from the Sx and S^ parts, respectively. Since there 
are no S^ or S^ terms in the Hamiltonian, we have tm,m±2 = tm,m±3 = 0. 

The recursion relation ( |2.5| ) may be interpreted as the Schrodinger equation of an elec- 
tron in a one-dimensional tight binding model. That is, we can consider the diagonal and 
off-diagonal terms as the on-site energy and hopping terms, respectively. Once this analogy 
is recognized, assuming tm,m±a{oi = 0, 1, or 4) vary slowly with m, we can treat the recur- 
sion relation within a continuum quasiclassical approximation or a phase integral method 



T3,IT3JT^. With this approximation we can define smooth functions 



taim) - ^itm,m+a + tm,m-a), = 0,1,4. (2.6) 



For the present model, t^'s are given by 



2/1 I \ _2\ 



to^m) = —m (1 + AiTTi 
Sh 



ti{m) = -^^S^-m^, (2.7) 



2 _2\2 



where we have used the approximation S{S + 1) ~ S"^. Introducing the DPI wavefunction 
within the semiclassical approximation 



Cm ~ . exp 

'f (m) 



i / q{m')dm! 



(2.8) 

we have the Hamilton -Jacobi equation 

E = Hsc{(l,fn) = to^rn) + 2ti(m) cosg + 2ti{m) cos4g, (2.9) 

and the transport equation 

, , dTisc 
v{m) = ^- 

= — 2sing(r?2)[ti(m) + 16t4(m) cosg(r7i) cos2g(m)]. (2-10) 



In Eqs. ( ^.81 ) and ( |2.10D , q{m) is a local, ??7,-dependent Bloch wave vector obtained by solving 



Eq. (|2.9|) for q for any given energy E. It is very useful to have a physical picture of these 
equations. For a given value of m, Eq. ( p.9|) gives an energy band E{q) which defines the 



classically allowed range of energies. In Fig. |l] we show possible E vs. q curves for our 
problem. At lower and upper edges of the band the transport equation shows that v{m) 
becomes zero because the slope dE{q)/dq is zero. This means the band edges are related to 
the classical turning points. These are not the only turning points, however. Such points 
are more generally defined by the condition that the velocity v{m) vanishes. This condition 
produces additional loci in E — m space, which we call critical curves, along with the m- 
dependent band edges. These curves are crucial to understanding how the oscillating tunnel 
splitting, i.e., the quenching effect, appears. 

B. Critical Curves 



L^ 



From Eq. ( p.lO[ ) the condition v{m) = is satisfied when q = 0, or q = it, ot q = q^,, 



where q^ is the solution of 

32t4(m) cos'^g=K(m) — 16t4(m) cosg=K(m) + ti(m) = 0. (2-11) 



Substituting these into Eq. ( p.9| ) we obtain the following energy curves for each of the three 
g's 



U-. 



'o{m) = toM + 2ti(m) + 2U{m), 

;(m) = toM - 2ti(m) + 2U{m), (2.12) 

U^{m) = tQ{m) + 2ti(m) cosg*(m) + 2ti{m) cos4g*(m), 

where Uo{m) = E{0,m), Un{rn) = E{n,m), and U^{m) = E{q^{m),m). Whenever a given 
energy E crosses one of these curves a turning point occurs. Various types of turning points 
depending on the characteristic of the critical curves have been analyzed in Ref . [|1^ . An 
interesting feature of this analysis is the existence of novel turning points inside the classically 
forbidden region, which is crucial for the quenching of spin tunneling. The recursion relation 
studied there was based on a spin Hamiltonian which includes terms up to second order, 
and there were only three critical curves to be considered. Here, we expect to have up to 
five curves, Uo{m), U-^im), and up to three U^{my?> from the cubic equation (|2.11|) . 



In order to proceed further, it is necessary to analyze the critical curve structure more 
closely, in particular, its dependence on h^. To do this, let us first compare Uo{m) with 
f/n-(m). From Eq. (2.7) it can be easily seen that UT,{m) > Uo{m) since tQ{m) < 0, ti{m) < 0, 
and t4(m) > for all \m\ < S. Thus, Un{m) can be the upper band edge. However, in order 
for this to be so we still need to prove that U-^im) > U^{m). This is not obvious. Indeed, 
since the Eq. ( p.ll| ) is a cubic in cosg,,,, it is possible to have complex solutions. These 
solutions will yield a complex U^:{m), which is not of interest because the Hamilton- Jacobi 
equation E = U^^ {m) can not then be satisfied. A careful consideration of the solutions of 
Eq. ( p.ll|) is therefore necessary. 



Defining x = cosg*, n = m/S, and using Eq. (2.7) for the t^'s, we can write Eq. ( p.ll| ) 
as 

fix) ^2x'-x- ^A^(l - M^)-^/^ = 0. (2.13) 



A sketch of the function f{x) is drawn in Fig. |[ This sketch incorporates the following 
easily verified properties of f{x): (i) /(O) < 0, (ii) /'(O) = -1, (iii) /'(±1) = 5 > 0, (iv) 
/(-I) < 0, (v) /(I) may be of either sign, (vi) f (±1/^/6) = 0, where f'{x) = df{x)/dx. It 
follows that a curve of type marked (a), characterized by one real zero of f{x) arises when 
hx is large, or when \m\ is large, and that a curve of type marked (6), characterized by three 
real zeros arises when h^ is small, or when |m| is small. Let us denote the largest zero by 
Xi, and the other two, when they are real, by X2 and x^ with X2 > x^. The corresponding 
values for q^:{rn) and U^^im) are denoted by g=Kj and Un{m), with i = 1, 2, or 3. It is obvious 
that xi > 0, and that — 1 < X3 < — 4g < X2 < 0. The first real root yields a positive value 
for cosg^Ki, but since we cannot say if xi is greater or lesser than 1, g=Ki may be real or pure 
imaginary. The other two real roots, when they exist, always yield real wavevectors g=K2 and 

q*3- 

The transition from one to three real roots occurs when f{x) has a double zero, i.e., f{x) 
and f'{x) both vanish simultaneously. It is easily shown that this condition is equivalent to 



hxc{m) = h 



h^ 




2\3/2 



(2.14) 
(2.15) 



The curve h^d^ and some special values of h^ are displayed in Fig. ^. The physical 
meanings of these values are listed in Table |I[ From the arguments of the previous paragraph, 
it follows that we will have three zeros when /i^. < hxc{tTi), and one zero when h^ > hxc{Tn). 
When hx < h^max, we can also ask for the points ±m,a{hx) at which we change from one to 
three real roots of f{x). These are directly given by solving Eq. (|2.14|) for hxdjn) = h^: 



rrir, 



S 



h,,. 



2/3' 



. '''i'max , 



1/2 



(2.16) 



Next, let us investigate whether U^:i{m) is inside or outside the classically allowed energy 
band. Since xi moves to larger positive values as \m\ increases (see Fig. |), we see that C/*i 
lies inside the band if xi < 1, i.e., \m\ < m*, where m* is such that /(I) = 0. Solving this 
equation we get 



m* 



S 



1- 



llxr 



2/3' 



1/2 



/27 



(2.17) 
(2.18) 



Clearly, nia < m^- 

Let us also explore whether the f/*j(r?7,)'s, when they are real, are larger or smaller than 
UT^{m) or UqIiti). We consider the following differences: 



UT,j{m) = UT,{m) — U^j{m) 

= 16t4(m)cosg^,j(l + cosg*j)^(-2 + Scosg^,^), 



(2.19) 
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= -16t4(m)(cos^ q^i - cos^ q^j) 3(cos^ g*i + cos^ q^j) - 1 , (2.20) 

Uoi{m) = Uo{m) - U^i{m) 

= 16t4(m) cosg=Kj(l — cosg=Kj)^(2 + 3cosg*j), (2.21) 

where i,j = 1,2, or 3 and we have used Eq. ( p.ll| ) to ehminate ti{m) in favor of t^lm) . 
From these equations, and using the facts that ti{m) > 0, plus [|24| 



11 1 , , 

—^ < cosq^s < -= < cosg*2 < 0, —= < cosg^i, (2.22) 

v2 v6 v2 



we find : 

1. When there is only one real root. 



f/,i(m) < f/o(m) < U^{m) (2.23) 



for all \m\ < S and /i^^max < hx < hxco- 
2. When there are three real roots. 



[/*i(m) < f/o(m) < f/*3(m) < f/*2(m) < U^{m) (2.24) 



for hxi < hxi 



U^i{m) < U^3{m) < Uo{m) < U^2{m) < U^{m) (2.25) 

for < h^ < hxi, where hxi is determined by Uo{m = 0, hxi) = U^:z{rn = 0, hxi), which from 
Eqs. (|2.21| ) and (p.22|) , is equivalent to cosg=K3(m = 0, hxi) = —2/3 . 



We can now list the various types of critical curve patterns that arise in our problem, 
and the corresponding ranges of the field hx- In the following, ?7_(m) and U^{m) denote 
the lower and upper bounds of the energy band, and Uf{m) and Ui{m) mean the forbidden 
and internal energies, respectively. 

In this case [/*2 and [/^.s are not real for any m, and q*i{rn) is imaginary, i.e., f/=Ki(m) is 
outside the band for all \m\ < S. The energy band E{q) is of the type in Fig. |l]a for all m, 
and the critical curves become 

U,i = Uf, f/o = f/-, U^ = U+ for \m\ < S, (2.26) 

which are shown in Fig. ^. 

L^aSe ii . flxmsix ^ '^x ^ '^xr' 

Now, f/*2 and f/*3 continue to be complex for all m, but g*i(m) is real in the central region 
\m\ < m=K. In this region, the energy band is as in Fig. |l|b, while in the outer region it is of 
the type in Fig. |l|a. Accordingly, the critical curves have the structure shown in Fig. ^ and 
can be written as 

f/*i = f/-, f/o = Ui, f/yr = f/+ for \m\ < m*, 

f/^i = Uf, f/o = f/-, U^ = f/+ for \m\ > m,. (2.27) 



Case III : < h^ < /i^max- 
There are now three m regions. In the outer region, \m\ > m^, [7^,2 and t/^,3 are still 
complex, U^i is outside the band, and E{q) has the shape in Fig. 0a. In the intermediate 
range m^ < \m\ < m^,, U^2 and U^^ continue to be complex, but U^i is inside the band, and 
E[q) has the shape in Fig. |l|b. In the central range, \m\ < rria, U^2 and f/*3 become real, and 
E[q) has the shape shown in Figs. |l|(c) (when h^i < h^ < h^max) and {d) (when h^ < hxi). 
The critical curves can be expressed as |2^ 

Uq, U^2, U^3 = Ui, f/*i = f/_, U^ = U+ for |m| < m^, 

Uq = Ui, U^i = U-, Uj, = U+ for TUa < \m\ < m*, (2.28) 

which are illustrated in Figs, ^a and jpD. 

As a matter of fact, we should distinguish two subcases in Case III. When h^ > hxi, as 
in Eq. ( p.24| ), the relevant critical curves are as in Fig. |]a. When h^ < hxi, as in Eq. ( p. 251 ), 
there is a range of m values in which t/*3 < Uq (see Fig. ^d). For the experimental parameters 
relevant to Mni2-ac, the field h^i is rather small, and the points rriQ, m^, rua, and rrii are 
all clustered tightly near m = S. This means that for the low lying states, there will be 
four turning points very close to one another, and the DPI analysis would have to be done 
using a quartic turning point formula, analogous to the quadratic turning point formula as 
discussed by Berry and Mount |^. Since we know the qualitative structure of the energy 
spectrum for fields as small as h^i, based on the arguments of Sec. Ill, e.g., this exercise is 
largely academic, and we have chosen not to perform it. This means that our analysis is not 
quite correct at very small fields, and this can be seen in Fig. |^ especially in the behavior 
of the splitting between the first excited pair of levels. As we shall discuss in Sec. Ill, this 
splitting is rigorously zero at h^ = 0, whereas we appear to find a zero at a slightly non-zero 
value of hx- 

As discussed in Ref. [|T^ the quenching of spin tunneling occurs when q{m) has a real 
part as well as an imaginary part inside the forbidden region. From the viewpoint of energy 
curves this happens when there is an energy curve inside the forbidden region. From the 
above analysis we can see that only f/*i(m) resides inside the forbidden region. For a given 
energy E such that f/omin ^ E < f/^imax, Q changes from pure imaginary to complex as 
m passes from the \m\ > rric region to the \m\ < rric region, where rric is the point where 
E intersects U^i{m), (for example see Fig. ^). When q becomes complex the semiclassical 
wavefunction in Eq. ( |2.8| ) oscillates with exponentially decaying or growing envelope. The 
quenching of spin tunneling arises from this oscillating nature of the wavefuction inside the 
forbidden region. 

We note here that for the experimental Mni2 parameters, the field /i^max is quite small 
(see the legend in Fig. ^, and so in the entire field range for Case III, even though there 
is a forbidden region turning point, the behavior of the ground state tunnel splitting is 
qualitatively similar to that for h^ = 0. The behavior of the splitting of the next two levels 
is more interesting, and as can be seen from Fig. |^, the DPI method does capture it, at least 
qualitatively, and perhaps even quantitatively. 



C. Tunnel Splittings 

We now calculate the energy splitting due to the spin tunneling between degenerate 
states in Mni2-acetate. In Ref. |I9|, tunnel splittings for five-term recursion relation have 
been obtained from Herring's formula. The final result is, however, quite general so that it 
can be applied to a recursion relation which includes more than five terms. Moreover, as 
we can notice from the above classifications, although the present nine-term case has more 
critical curves the possible types of the turning points are all included in those discussed 



in Ref. [113] , and no new type of turning point emerges here. Thus, we can directly apply 



the formula for the tunnel splittings obtained in Ref. |]T9| to the present problem. Since our 
calculation is based on this formula we quote the main results here. The tunnel splitting for 
nth pair of states is given by 



nl \l TV 



(2.29) 



where 



An 

F 

Qi 



mo 
rric / X 

Xo + (^ + 2)^0X0 ) dm, 
2Muo (mo — rric) 



ririLLlll I . 



(2.30) 



X exp ( -2 Qi + ujq 



-mo 




dm. 



dK{mi, e) 



Here, k and x ^'^'s the imaginary and real parts of complex q, respectively, and 

Kq = K{m, e = 0); k'q - 

Xo = x{m,e = 0); x'o - 
Bq = cos g(m, e = 0);Bq 



de 
dxim, e] 



e=0 



^e .=0 

d COS q{m,, e) 



f2.311 



de 



e=0 



with e = E — f/_(— mo). In these equations, ±m,c are not quite the turning points of the 
previous subsection, in that they are not the point where [/^,i(m) equals the true energy 
En of the nth pair of levels. Rather, they are the points where f/=Ki(m) = f/_(±mo), which 



corresponds to setting E = f/_.(mo), i.e., e = 0. The reason is that the formula (|2.29|) 
incorporates expansions of various phase integrals in the energy difference 



En - U-{mo) = {n + -)uJo, 



(2.32) 



which is of order (1/5) compared to the energy barrier, as long as n ^ 5. This is why m^ is 
modified, and also why the primary phase integral for the Gamow factor To runs from — m-o 



9 



to rrio, the ininiina of UQ{m), rather than between the points where Uoim) = En- Since all 
energy curves are a function of both m and h^, these points still depend on h^, which in 
turn makes the A„ depend on hx- 

The mass M and frequency ujq in Eq. (|2.3CI|) are obtained by approximating ?7_ {m) near 
its minima by a parabola, i.e., U-{m) = E + ^mujQ{m ± mo)^. For m = —mo we find. 



M. ' 



, 5 



2ti(— mo) + 32t4(— mo) 
cJq = -2[ti(-mo) + l6U{-mo)] 



dm? 



(2.33) 

ni=— mo 



The application of the formulas ( |2.29|) -( p.33| ) cannot be carried out in closed form all 



the way, and we must resort to numerical methods. We explain the principal steps in our 
numerical calculation below. 

In step 1, we must find ±?7io, and [/_(±mo)- For our problem we discover that U-{m) 
is always given by Uo{m,) near the classically allowed minima. The equation for the minima 
can be reduced to another cubic, 

S^hl = 4(^2 _ ^)[2(Ai - 2\2)y + 1 + 4X28^]', (2.34) 

where y = m?. For th parameters Ai and A2 of interest to Mni2, and h^ < h^co, all three 
roots of this cubic equation are real, but only one is positive. This root gives us mo, and 
substitution of this value into Eq. (|2.12| ) for t/o(m) gives E, and Eqs. ( p. 331) then give M 
and uq. 

Step 2 is to obtain the points ±m,c given by the roots of the equation 

f/,i(m) = t/_(mo). (2.35) 



As discussed after Eq. (|2.32| ), up to terms of relative order (1/5*), the points ±m,c are 
the actual turning points for the low lying energies. Note that it is f/=Ki which appears in 
Eq. ( |2.35D since this is the critical curve that lies in the classically forbidden region. 



To solve Eq. (|2.35|) numerically, we first solve Eq. (|2.11|) for the function cosg*i(m), which 
can be done in closed form. This solution is then substituted in Eq. ( ^.12| ) to obtain [/*i(?7i). 
The entire procedure can be implicitly implemented in the numerical routine. The same 
holds for dU^i{m) / dm. Since U-{mQ) is known from step 1, any of the standard root-finding 
methods — Newton-Raphson, bisection, secant etc. — can be applied to Eq. ( |2.35| ). 

Step 3 is to find q{m,), in particular its real and imaginary parts K,o{m,) and xo{^- This 



is done by solving the Hamilton- Jacobi equation ( |2.9| ) with the energy E found in the first 
step. The problem amounts to solving a quartic equation in cos q and making sure that one 
has the correct solution, which can be done easily by making use of the properties that we 
have found above. Thus in the region rric < |m| < mo, there are two solutions of the form 
in (with K, real), and two of the form ir — in. We discard the latter, and of the former select 
that one which continuously tends to as ?7i ^ ±mo. In the region |m| < rric, the solutions 
can be written as m ± x, and as vr — (m ± %), where % ^ as m — > ±mc- We discard the 
latter two, and read off k,{= kq), and x{= Xo) from the imaginary and real parts of the first 
two. Note that both kq and Xo are taken to be positive. 

10 



Step 4 is to find the e partial derivatives Kq and x'o, in effect dq{e, m)/de. {B'q is directly 
obtainable from Kq.) We differentiate the Hamilton- Jacobi equation with respect to E: 

-2[tising + 4t4sin4g]-^ = 1. (2.36) 

Since q{m) is found in step 3, this equation gives dq/de for any m directly. 

We now have all the ingredients needed to evaluate the one-dimensional integrals Fq, A„, 
F, and Qi. This is a straightforward numerical procedure. The only point worth noting is 
that the integrand for Qi is non-singular at m = — ttiq, and behaves, in fact, as (m + mo). 

In Fig. 1^ we show the tunnel splittings for first three pairs of states as a function of the 
field parameter h^. For comparison we also plotted the tunnel splittings obtained from exact 
diagonalization of the Hamiltonian. 

lYionv the results we observe several interesting features. First, as anticipated the tunnel 
splittings are completely suppressed at certain values of hx- The overall pattern of zeros, 
their number, and the dependence of this number on n, the pair index, is understandable 
on general grounds as we shall discuss. What is surprising is how regularly spaced these 
values of hx are. For the first pair of splittings, e.g., the intervals between successive zeros 
decrease by 2 or 3% only, and the last interval is 92% of the first. For the next pair, A2, the 
last interval is 95% of the first. The mean interval between zeros for the first three pairs is 
A/7o - 0.93T, Aifi ~ 0.85T, Ai^o - 0.79T. 

The regularity of the zeros means that the phase integral A^ decreases almost linearly 
with hx- (From Fig. |^, Gamow factor Fq also appears to be quite linear in hx.) While this 
variation is clearly expected to be smooth, we have no a priori way to judge how linear 
it will be. A similarly strong regularity of quenching intervals is experimentally discussed 
in Feg. The simplest model Hamiltonian for Fes entails only second order terms in the 
components of the spin operator, and in this model, the spacing of zeros is exactly equal 
28| , but to describe actual Fcg, one must add fourth-order terms. These terms change the 



spacing significantly, but still seem to preserve its regularity. It would be interesting to find 
a physical argument for this feature, which appears to be somewhat general. 

Second, aX hx = Q the tunnel splitting alternates between zero and nonzero as the level 
number goes up. This is due to the fourth-order terms in the spin Hamiltonian. This term 
causes the tunneling in the Mni2 system, but it also restricts the transitions to the case that 
the difference. Am = |m — m'|, between levels E^ and E^n' is a multiple of 4. For tunneling 
between degenerate states E^ and E^m this requires the condition 2m = 4p, where p is an 
integer. Therefore, there is no tunneling between +m and —m when m is an odd number. 

To give a more detailed argument of this point, we note that when Hx = Hz = 0, because 
of the S^ terms, the Hamiltonian can be divided into the following subspaces for S* = 10. 

Vi = {-10, -6, -2, +2, +6, +10}, 

\/2 = {-9,-5,-l,+3,+7}, 

\/3 = {-8,-4,0,+4,+8}, (2.37) 

V4 = {-7,-3,+l,+5,+9}, 

where the numbers in brackets give the m quantum numbers. The subspace Vi contains 6 
levels, which form 3 pairs split by tunneling due to the CS^ terms. The space V3 contains 
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5 levels, of which ±8, and ±4 are split by tunneling, and m = is isolated. There is no 
degeneracy amongst the 5 levels in space V2, but because of time reversal, this space is 
isomorphic to V4, and we therefore conclude that in the full spectrum of 7i, there should be 
five pairs of strictly degenerate levels, corresponding approximately to m = ±(2n + 1) with 
integer n. 

Third, there are five quenching points in the ground and first excited states tunnelings, 
and the number of quenching points decreases as the level number goes up. For the ground 
state tunneling the allowed number of quenching points can be explained qualitatively by 
thinking in terms of instantons. Since the instanton approach also gives a good geometrical 
structure to the quenching we present it in a separate section. 

III. QUALITATIVE TREAMENT 

In this section we give a qualitative treatment of the spin tunneling in Mni2 using in- 
stanton methods. Adding a hard axis field to Eq. ( |1.2| ), the Hamiltonian becomes 

H = -AS^, - BSt + C{SX + St) - g^iBH^S,. (3.1) 

The instanton method is based on spin coherent state path integrals. In the spin coherent 
state representation the anisotropy energy corresponding to the Hamiltonian is given by 

nc{a,(3) = {n\n\h) 

= —AS'^ sin^ a sin^ /5 — BS'^ sin^ a sin^ (3 

+2CS"*(cos'^ a + sin^ a cos^ /3 — 6 sin^ a cos^ a cos^ P) — qjibHtS cos a, (3.2) 

where a and (3 are the polar and azimuthal angles with respect to x, i.e., 

cos a = n ■ X, 
sinacos/5 = n ■ y (3.3) 

sin a sin /5 = n ■ z. 

The energy (p.2|) exhibits two degenerate minima at fij = (a,/?) = (Q;o,vr/2) and fi/ = 
(a,/3) = (ao;~7r/2), where Oq = pV^ for H^ = 0, and decreases smoothly to as H^ is 
increased. The level splitting due to tunneling between these minima can be obtained from 
the imaginary time propagator 

Kfi = {hf\exp[-nT]\hi) 

= Jv[h]exp[-SE[n{r)]], (3.4) 



where 



with 



SE[n{T)] = -tSA\h{T)]+ f nc{a,(3)dT, (3.5) 

Jo 

A[fi(T)]= [ (1 -COS a) p(T)dT, (3.6) 

Jo 
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is the Euclidean action and complex in general. Here, the boundary conditions are n(0) = fij, 
n(T) = hf. Geometrically, for a closed path, the integral in Eq. ( p.6|) can be interpreted 
as a surface area on the unit sphere enclosed by this path, which can be verified by Stokes 
theorem. In the large spin limit the path integral can be approximated by the sum of 
all contributions from paths that minimize the action, that is, the instanton paths. The 
instantons for the present model are not simple because of the fourth-order terms. However, 
we can construct a qualitative argument to find the quenching effect without performing 
explicit calculations . Since the Euclidean action Se has both real and imaginary parts we 
can generally express the ground state tunnel splitting as 

A = Y,Dje~^'''e'^'^, (3.7) 

j 

where j labels the various instantons, Srj, Sjj are the real and imaginary parts of the 
instanton action, respectively, and Dj are prefactors. With these ingredients we now discuss 
how the quenching appears in the present model. 

Let us first consider the case when H^ = 0. Since the energy has four- fold symmetry, 
an argument of von Delft and Henley can be applied ||^. If n(T) is an instanton path, so 
is 7^z(7r/2)n(r), where 7^z(7r/2) is a rotation through 7r/2 about z. Keeping in mind that 
n(r) is complex, when we project onto the real unit sphere, there are four saddle point 
paths passing through each of the four medium directions. Because of symmetry, each has 
the same real contribution to the action integral Sr. However, since their azimuths about 
the easy axis are different, the imaginary part of the action, i.e., the phase 5*/, will not be 
the same. From the geometrical meaning of the integral in Eq. ( p.6|) , the phase difference 
between two instanton paths equals S times the surface area on the unit sphere enclosed by 
these instanton paths. 

To visualize the interference effect we map the two-sphere onto a plane, as in an ordinary 
atlas (see Fig. ^). The hard axes are mapped onto four equally separated points lying on the 
equator, and the points exactly halfway between these correspond to the medium axes. Thus, 
the real projections of the instanton paths can be drawn as curves which start from +z, pass 
through the medium points, and end at —z. The area enclosed by two adjacent instanton 
paths equals vr, since the sphere is equally divided into four paths by the instantons. Thus, 
the phase difference between adjacent paths becomes Sir. Choosing the phase of path 1 
as the base, we can perform the summation in Eq. ( p.7|) . Recalling that by symmetry the 
contribution from real parts of the instanton actions are all same, as are the prefactors Dj, 
we have 

= ADe-^^e'^ cos(7rS) cos {^\ , (3.8) 

where 7 is an irrelevant phase. This result gives us two quenching conditions. From the 
factor costtS", we obtain the quenching of spin tunneling for half-integer S, which is just the 
Kramers degeneracy effect. The second cosine implies that the ground state spin tunneling 
is quenched for odd integer spins, i.e., S = 1,3,5, etc, and so A is non-zero only for S = 2p, 
where p is an integer [^. 

We now consider the case with H^ ^ 0. Since the field is assumed to be applied along 
the -|-x axis, both easy and all four medium axes move close to +x axis. Thus, the two 
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dimensional picture becomes one as shown in Fig. ^ The four-fold symmetry is now broken, 
but there are two pairs of instanton paths surrounding the +x axis: {a, a') and (b,b'). The 
real parts of the instanton actions in a pair are same, but different between the pairs. The 
phase differences in each pair are the areas enclosed by each pair of instanton paths (the 
small and large oval regions in Fig. ^ and are dependent on the field Hx. If we choose 
the straight hne joining +z' to —z' as a reference, Sja' = —Sja, Sjb' = —Sib, so that the 
summation in Eq. ( |3.7| ) can be performed as 



A = DaC-^'"' \e-'^'^ + e-'^i-'] + D^e 



-Se 



^~iSjt _j_ ^—iSjy 



= 2D,e-^- cos ^MIA + 2D,e-^- cos ^^^M^, (3.9) 

where SRa, Sub are the real parts of the instanton actions in each pair, and Aa{Hx), Ab{Hx) 
are the areas enclosed by the pairs {a, a') and {b,b'), respectively. For H^ > the saddle 
points through which the paths (a, a') pass are lower than those for (6, 6'), which means that 
SRa < S^b- The main contribution to A in Eq. ( p.9| ) then comes from the first term, and 
we can neglect the second term. The quenching of the ground state tunnel splitting thus 
arises when Aa{Hx) ~ (2n + l)7r/2S', where n is a non-negative integer. To see how many 
quenching points are allowed we note that Aai^H^) < ^^(O), where Aa{0) = n (the area 
enclosed by the two paths 1 and 4 in Fig. ^. From this condition we find n < {S — l)/2. 
For S* = 10 there are thus five values of H^ at which the quenching appears. 

The instanton argument provides another way of seeing that the region of very small 
(but non-zero) H^ is special. Exactly at H^ = 0, four instantons are important, but for 
large H^ only two are important. There must therefore be a regime of small H^ where we 
make a smooth transition between these two behaviors. The width of this regime can be 
quite small since Sjib > Sj^a as soon as H^ j^ 0, and these actions appear in the exponents in 



Eq. (|3.9|), so that the difference {Sub — S^a) is amplified. We thus have another way of seeing 
why the formula (|2.29|) fails near h^ = 0. It contains only cosine factor, and is effectively 
ignoring the second term in Eq. ( p.9| ). 

IV. SUMMARY 

In this paper, we have used the DPI method to study tunneling in Mni2, especially its 
behavior with a hard-axis field, which is expected to show oscillation as in Feg. The recursion 
relation now has nine terms, complicating the analysis. There may be up to five critical 
curves, which leads to many more turning points. The DPI method still works, however, 
even though the phase integrals and integrands must be evaluated numerically. But, the 
numerical procedures required are simple, and involve only root finding and integration in 
one variable. Except for some special narrow field regions, where two or more turning points 
merge, the DPI analysis based on linear turning point formulas is extremely good, and agrees 
with exact numerical results quantitatively. 
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FIG. 1. The different types of energy bands possible for recursion relation ( p.5| ) with wavefunc- 
tion ( |2.8| ). Note UT^{m) = [/+ in all cases, (a) When hxr < hx < hxco with |m| < S, and when 
hxc{^ < hx < hxr with |?TT,| > rn,*. In this case, UQ{m) = C/_, and [/*i(r?T,) does not appear since q^i 
is imaginary. (6) When hxci^n) < hx < hxr with nia < \m\ < m^,. Here, C/=Ki(m) = [/_, f/o("^) = Ui. 
(c) When /i^^j < hx < ^xmax with |m| < nia, and when < /i^; < /i^.^ with rrii < \7n\ < nia- (d) 
when < hx < hxi with |?tt,| < rtii. Note that, in both (c) and (d), C/=Ki(m) = C/_, and Uo{m), 
f^*2("^)) and U^3{m) are inside the band and thus all denoted Ui. 

FIG. 2. Sketch of the cubic function /(x) for (a) large /i^, or large |r?T,|, (6) small hx, or small 
|?TT.|. Note that, there is one root in (a), but three roots in (6). The transition from type (a) to type 
(5) occurs when f{—Xm) = f'{—Xm) = {xm = 1/V6), which gives the curve hxc{rn) in Eq. (2.14). 



FIG. 3. The curve hxd'm) and some physically meaningful values of /i^'s. In the inset we list 
these values computed with experimental numbers for Ai and A2 for Mni2-ac. Points at which hx 
intersects the curve hxciin) are m = ima- 

FIG. 4. The critical curves for Case I. At points ±m-o, Uq has minima, and the points inic 
denote the intersection between E and [/*i. Note U^i = Uf and Uq = U- for all \m\ < S. For a 
given value of E, q becomes complex for |m| < rUc which lies inside the classically forbidden region. 
In this region the semiclassical wavefunction Cm oscillates with decaying or growing envelope. 

FIG. 5. The critical curves for Case II. ±771^, are the points where UQ{m) = U^i{m) [and 
dUQ{m)/dm = dU^,i{rn)/dm]. U^,i is the lower band edge in the central region \m\ < m^ and 
forbidden in the outer region |m| > m^,. 

FIG. 6. The critical curves for Case III. (a) When hxi < hx < ^xmax and (b) when < hx < hxi- 
There are five critical curves. Note, however, that C/*2 and U^:^ appear only in the region \m\ < ma 
because they are complex outside this region. 

FIG. 7. Tunnel splittings A„ between first three pairs of levels as a function of the field pa- 
rameter hx . The dotted and solid curves are obtained from exact numerical diagonalization of the 
Hamiltonian and the DPI method, respectively. 

FIG. 8. Two dimensional picture of instanton paths when Hx = 0. The points zizx, ity are the 
hard axes, and ttzj's represent the medium axes. Dotted lines denote the real projections of the 
instanton paths. 

FIG. 9. Two dimensional picture of instanton paths when Hx 7^ 0. The points itz' represent 
the new easy axes. The instanton paths are again denoted by dotted lines. Note that the areas 
enclosed by each pair of instanton paths are shrunk due to the field. 
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TABLE I. Physical meanings of the special hx's 



Coercive field above which no tunneling exists. 



The value below which the wavefunction can have 
an oscillating part inside the forbidden region. 



The value above which g^<i becomes real. 



The maximum value of the curve hrrim) 



The value at which [7*3(0) intersects Uo{0). 
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